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Abstract: Recent tentative experimental indications, and the subsequent theoretical 
speculations, regarding possible violations of Lorentz invariance have attracted a vast 
amount of attention. An important technical issue that considerably complicates de- 
tailed calculations in any such scenario, is that once one violates Lorentz invariance the 
analysis of thresholds in both scattering and decay processes becomes extremely subtle, 
with many new and naively unexpected effects. In the current article we develop several 
extremely general threshold theorems that depend only on the existence of some energy 
momentum relation E(p), eschewing even assumptions of isotropy or monotonicity. We 
shall argue that there are physically interesting situations where such a level of gener- 
ality is called for, and that existing (partial) results in the literature make unnecessary 
technical assumptions. Even in this most general of settings, we show that at thresh- 
old all final state particles move with the same 3-velocity, while initial state particles 
must have 3- velocities parallel/anti-parallel to the final state particles. In contrast the 
various 3-momenta can behave in a complicated and counter-intuitive manner. 
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1 Introduction 



The OPERA collaboration's announcement of tentative but statistically significant in- 
dications of "faster than light" neutrinos [1] (see also earlier more tentative results from 
the MINOS collaboration [2]) has ignited a firestorm of theoretical speculation. Over 
150 theoretical articles have been generated in some 9, weeks. Notable contributions 
include [3-21]. In addition to these very recent efforts, it is important to recognize that 
there is an older and extensive literature placing significant experimental and observa- 
tional bounds on any possible violation of Lorentz invariance. See for instance work 
by Coleman and Glashow [22, 23], Jacobson and collaborators [24-33], and especially 
the Living Review by Mattingly [34]. Other theoretical frameworks for characterizing 
possible violations of Lorentz invariance include those of Nielsen and collaborators [35- 
39], Kostelecky and collaborators [40-48], and the flat-space non-gravity framework 
developed by Anselmi [49-61]. Additionally, the Horava gravity framework [62] natu- 
rally includes Lorentz violation [63-68] , though quantitatively the relationship between 
Hofava gravity and the OPERA results seems somewhat strained [14]. In this arti- 
cle we shall not directly address the phenomenology of the OPERA-MINOS results. 
Instead we shall prove some very general theorems on Lorentz violating processes, the- 
orems that should be borne in mind whenever one is attempting to step outside the 
framework of standard special and general relativity. 

One of the key results in the literature devoted to possible violations of Lorentz 
invariance is that the normal intuition one develops regarding threshold phenom- 
ena requires significant modification. See for instance the articles by Coleman and 
Glashow [22, 23], and Jacobson, Liberati, and Mattingly [26, 27] - - and the more re- 
cent follow-ups by Cohen and Glashow [5], and Liberati, Mattingly, and Maccione [19], 
focussing specifically on the OPERA results. In this article we will generalize the 
analysis of threshold phenomena presented in those articles. We shall consider both 
single-particle decay processes, and two-particle scattering processes (possibly inelas- 
tic), taking care to make an absolute minimum of technical assumptions — thus greatly 
generalizing previous analyses. Specifically: 

• We will explicitly assume a normal spacetime manifold based on IR 4 . 

(This excludes, for instance, both non-commutative spacetimes and certain ver- 
sions of DSR [69, 70]. This is a purely pragmatic decision based on the fact that 
we want to be able to say something reasonably concrete.) 

• We shall explicitly assume conservation of both energy and momentum. 

(This is again a purely pragmatic decision based on the fact that we want to be 
able to say something reasonably concrete.) 
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• Furthermore we shall explicitly assume a Hamiltonian/Lagrangian framework, so 
that in view of Noether's theorem (combined with energy-momentum conserva- 
tion) we are working in a homogeneous spacetime. Specifically, we assume the 
free-particle energy to be some function of the 3-momentum, E(p), and that this 
can be related to a 3- velocity via Hamilton's equations 



(This is again a purely pragmatic decision based on the fact that we want to 
have a sufficiently well-defined framework in which to be able to say something 
reasonably concrete.) 

• We shall eschew any particular functional form for E(p), though we will impose 
smoothness and differentiability constraints as needed. 

• We shall explicitly not assume isotropy. 

— Even if physics happens to be isotropic in the preferred ("aether") frame 
implicit in many specific Lorentz violating theories, there is no particular 
reason to assume isotropy of the energy- momentum relation in generic in- 
ertial frames. And since we do not necessarily know what the observer's 
3-velocity is with respect to the preferred frame, it is more useful to develop 
threshold analysis for generic observers in an explicitly observer-dependent 
manner. 

— Even as early as the 1980's attempts were made to take lattice physics 
seriously as a physical cutoff — with proton decay taking place with outgoing 
decay products preferentially aligned along the principal axes of the universe. 
In such a situation one would not have isotropy even in the preferred frame. 

— Many of the "analogue spacetime" models permit energy-momentum rela- 
tions that have odd and possibly anisotropic behaviors at ultra-high ener- 



— Consider a generic quasiparticle propagating in a generic atomic lattice. 
(For instance, a conduction-band dressed electron.) The band structure will 
typically not be isotropic, even in the rest frame of the lattice. So if you are 
a condensed matter physicist, you will have no choice, you will simply have 
to acknowledge that non-isotropy of the energy-momentum relation is quite 
common, which will unavoidably influence your ability to analyze reaction 
thresholds. 



dE 




dp ' 



P 



0. 



gies [71, 72]. 
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• We shall also be extremely cautious concerning "monotonicity" assumptions - 
carefully formulating an appropriate concept of monotonicity, and carefully ana- 
lyzing what can and cannot be extracted from such an assumption. 

Because of the generality of these assumptions, our results will have considerably wider 
validity than the results currently extant in the literature. Even in this most general 
of frameworks, several rigorous theorems can be extracted. We shall show that at 
threshold all final state particles move with the same 3-velocity, while initial state 
particles must have 3-velocities parallel/anti-parallel to the final state particles. In 
contrast the various 3-momenta can behave in quite complicated and counter-intuitive 
fashion, and the 3-momenta need not even be collinear. 

2 General background 

In the Coleman-Glashow analysis [23], most of the discussion is explicitly limited to 
the rather special case of single-particle decay processes where the initial and final state 
particles (i G {1,2,... , n}) all have energies of the form 

Ei( Pi ) = ^So + IIPilfc?. (2-1) 

Here the "speed of light" can be particle dependent. (See especially equation (2.19) in 
reference [23]. A similar assumption is implicitly made in [22].) But it is reasonably 
clear that much of the discussion of thresholds in [23] would work for any generic 
Ei(pi). In contrast, in the OPERA-related analysis of [5], this specific choice of energy- 
momentum relation is implicit, not explicit, but is absolutely essential to that discussion 
— see [19] for a generalization. 

The Mattingly, Jacobson, and Liberati threshold analysis [27] focusses on energy- 
momentum relations that are (even in their most general setting) taken to be both 
isotropic and monotonic. Specific examples are taken to be of the form 

Ei{Pi) = y^HWP^i {\\pi\\c) n El- n . (2.2) 

Here rji is a dimensionless parameter and E* is an energy scale characterizing the devi- 
ations from Lorentz invariance. The special case n = 2 corresponds to the Coleman- 
Glashow energy-momentum relation. If rji < and n > 2 then these energy-momentum 
relations can in principle exhibit a maximum — the energy "saturates" — this is a spe- 
cific example of a much more general phenomenon: 
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In lattice QFT regularizations energy-momentum relations are typically of the 
form 

*m = Mt)V(^) + ^(T) + ^(x)}- ^ 

or, (for "massless" particles), 



and typically exhibit a maximum energy E max = \/ Eq + 3h 2 c 2 /a 2 in terms of 
the lattice spacing a. (One usually considers "small" momenta, ||p||a//i <C 1, 
where Lorentz invariance is approximately recovered. Herein the focus will be on 
deviations from Lorentz invariance.) 

Qualitatively similar effects occur for quasiparticles propagating through atomic 
lattices — momentum space [the first Brillouin zone] is now compact, and so 
(assuming continuity) the energy will be bounded by some maximum. 

In DSR-inspired models, insofar as they can be incorporated into the current 
framework, one often has individual particle energies saturating at or around the 
Planck energy [69, 70]. 

For a specific example (physically unmotivated but mathematically tractable) of 
saturation behaviour one might take 



E(p) = J El + El tanh ( W!^) . (2 . 5 ) 



At low momentum this is approximately Lorentz invariant 



£(p)« ^E 2 + \\p\\ 2 c 2 + 0(\\p\n, (2.6) 
but at high momentum it exponentially saturates 



E(p) « + El + 0(exp(-2\\p\\ 2 cyEl)). (2.7) 

Neither the usual Lorentz invariant energy-momentum relation, nor even the 
Galilean invariant energy-momentum relation, saturate with a maximum energy. 
For that matter, neither does the Coleman-Glashow energy-momentum relation 
(2.1) saturate. 
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Since the whole point of this article is to generalize threshold analysis as far as possible 
to situations where Lorentz invariance is broken, we will for generality entertain the 
possibility of energy-momentum relations that saturate to some maximum energy - 
with the understanding that this maximum energy might, in specific situations, be 
infinite. With this general framework in place we are now ready to begin detailed 
analysis. 



3 Cautionary comments 

Some cautionary comments are in order: 

• For the lattice-like energy-momentum relation of equation (2.3) we have 
BE he ( . (2p x a\ . (2p y a\ . (2p z a\\ 



dp 2Ea \ \ h J' \ h J \ h 

The key point is that 3-velocity v and 3-momentum p need not be parallel. 
Additionally v can exhibit non-trivial zeros for non-zero momentum p, and even 
once one specifies a particular particle the inverse function p(v) can easily be (and 
typically is) multivalued. Such phenomena are not limited to the specific energy- 
momentum relation of equation (2.3), but rather are generic to any quasiparticle 
propagating through a regular lattice (for example, a conduction-band dressed 
electron) . 

Non uniqueness of the inverse function p(v) is also generic for (higher than 
quadratic) polynomial or rational polynomial energy-momentum relations — it 
is the unique invertability of the Lorentz invariant energy-momentum relation 
E = y/E% + ||p|| 2 c 2 that is non generic in this regard. (Details depend on the 
precise values of the coefficients as the potential for multi-valued behaviour de- 
pends on the root structure.) 

For the tanh-like energy-momentum relation of equation (2.5) we have 

v = ^=ese^(^)^. (3.2) 



dp V E * ) E 

While 3-velocity v and 3-momentum p are now parallel, zero 3-velocity can corre- 
spond either to zero 3-momentum or to infinite 3-momentum (with finite energy 
\J El + E^). Low- velocity physics can thus be grossly misleading — two particles 
with the same 3-velocity may have wildly differing 3-momenta. Such phenomena 
are not limited to the specific energy-momentum relation of equation (2.5), but 
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rather are generic to any situation where the energy saturates as a function of 
3-momentum. Note in particular that the energy-momentum relation of equation 
(2.5) is monotonic — monotonicity is not enough to prevent this sort of behaviour. 
Similar behaviour also occurs whenever lirni|p|i_ >0O dE/dp — > 0, corresponding to 
a sub-linear asymptotic growth in the energy-momentum relation. 

Other unusual possibilities include energy minima occurring at non-zero 3-momen- 
tum, (by definition an energy minimum must always occur at zero 3- velocity). 
Let n be an arbitrary unit vector and consider for instance 



E = ^JEl + | \p\ | 2 c 2 + k 4 (p ■ n) 4 + k 6 (p ■ n) 6 . (3.3) 

This energy-momentum relation is not only anisotropic, but by taking k 4 < and 
k 6 > one can arrange for a global minimum energy at some p min = p m i n n ^ 0. 
The 3-velocity is 

_ PC 2 + 2/c 4 (p ■ n) 3 n + 3fc 6 (p ■ n) 5 n 

and p and v are generally not collinear (unless one happens to be considering 
motion parallel or perpendicular to the preferred axis n.) 

As a final pedagogical example consider the isotropic energy-momentum relation 



E = sJEl + ||p|| 2 c 2 + h\\p\\* + fc 6 ||p|| 6 (3.5) 

By taking k± < and ke > one can arrange for a global minimum energy at 
some finite ||p m in|| = Pmm ^ 0. The 3-velocity is 

||p||c 2 + 2fc 4 ||p|| 3 + 3fc 6 ||p|| 5 A 
v = p. (3.6) 

In this situation p and v are generally collinear, but whenever there is a global 
energy minimum at some finite ||p m in|| = Pmin 7^ there will also be a non-empty 
range of momenta for which p and v are anti-parallel. 

These are merely five specific examples of the unusual behaviour one might potentially 
encounter, and the types of issues we shall potentially need to consider in our analysis. 

4 Decay thresholds 

Consider the decay process 

X ^X 1 + X 2 + ...X n , (4.1) 
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where for each individual particle we have the 4-momenta 



Pi = (Ei(p) iPi ). 



(4.2) 



We shall now study the kinematics of this decay process. 



4.1 Kinematically allowed region 



Let us define 



n 



£out(Pl,f>2,- • • ,Pn) = ^Eiipi) 



(4.3) 



i=l 



and 



^out(Po) = \ (Pl,P2, ■ ■ ■ ,Pn) ■ ^Pi = PO 



(4.4) 



Then P ut(Po) is a collection of three individually connected hyperplanes in M. n , one 
hyperplane for each Cartesian component of po, corresponding to the set of all pos- 
sible outgoing 3-momenta for fixed total 3-momentum p . Thus V out (po) is a 3n — 3 
dimensional plane (affine subspace) of co-dimension 3 in M 3n , and is both convex and 
connected as a subset of M. 3n . But the individual E^pi) are by assumption differentiable 
and continuous, so E out (p 1 ,p 2 , . . . ,p n ) is also differentiable and continuous. In partic- 
ular, since "P ou t(Po) is connected, this implies the image E ont (V ut(Po)) is a connected 
interval in M. 

For specified initial 3-momentum p , the decay process (4.1) is kinematically al- 
lowed if and only if 



That is, the decay is allowed if and only if among the set of all possible output 3- 
momenta {pi,P2, ■ ■ ■ ,Pn} that conserve total 3-momentum, there is at least one con- 
figuration that also conserves total energy. We could also phrase the kinematically 
allowed region in terms of an allowable set of output momenta by considering the 
inverse image 



a set which, for given po, may or may not be empty. Alternatively one can ask the 
question 




(4.5) 




(4.6) 




(4.7) 



But these approaches in terms of inverse images can be somewhat clumsy. 
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As a more practical way to better characterize the kinematically allowed region, it 
is useful to introduce the two quantities 

E min (p ) = min E out (V OVLt (Po)) = min <| ^Ejjpj) 

[ i=i 

and 

E max (p ) = max E out (V out (p )) = max <^ y^(pQ 

I i=i 

A more technically precise statement would use the concepts of supremum and infimum, 
but as long as we understand that statements made below might sometimes have to 
be interpreted in terms of suitable limits, such a level of precision is, for our purposes, 
unnecessary. 

Then the decay process (4.1) is kinematically allowed if and only if 

£min(Po) < £ (Po) < ^max(Po)- (4.10) 

A decay threshold is now defined to be the edge of the kinematically allowed region. 
Specifically, an enabling threshold is defined by the condition 

E min (p Q ) = E (p ), (4.11) 

and a saturation threshold is defined by the condition 

£ (P0) = ^max(Po)- (4.12) 

These thresholds are typically 2-surfaces in 3-momentum space. 

• Note that -E m ax(Po) might trivially be infinite if any one of the energy- momentum 
relations does not saturate at large 3-momentum, in which case no useful upper 
bound, and hence no saturation threshold, would be obtained. This is the case 
for instance in standard special relativity, in standard Galilean kinematics, in 
the Coleman-Glashow energy-momentum relation (2.1), and for rj > in the 
Mattingly-Jacobson-Liberati energy-momentum relation (2.2.) 

• Note that in standard special relativity the enabling threshold is also trivial - 
one need merely go into the centre-of-momentum frame to see that the decay is 
kinematically allowed if and only if 

n 

m > ^m,. (4-13) 

n=l 

That is, decay thresholds are trivial in the case of exact Lorentz invariance, (see figure 1) 
and only become interesting if there are deviations from Lorentz invariance. 



:J>i = PoL (4.8) 



»=i 



Po 



(4.9) 



i=l 
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Figure 1. The kinematically accessible region for a Lorentz invariant particle of mass m, 
and a Lorentz invariant particle of mass 4m, decaying to two identical particles of mass m. 
Note absence of decay thresholds: The process is either allowed or forbidden in a momentum- 
independent manner. 



4.2 Thresholds in momentum space 

For a graphical understanding of the situation it is useful to pick some (arbitrary but 
fixed) direction p in momentum space, write p = p p, and for each direction p 
consider the three curves: 




(4.14) 
(4.15) 
(4.16) 





Note that if the individual energy-momentum relations are isotropic, (rotationally in- 
variant, spherically symmetric), then these curves C+/o/- will be independent of the 



direction p. If Lorentz invariance is violated, isotropy would at best occur only in the 
preferred (aether) frame, so in general it is safer to not make any such assumption. 

The kinematically accessible region for the decay products (assuming only conser- 
vation of 3-momentum) is the region between the curves C_ and C+. Kinematically 
allowed decays correspond to that portion of Cq that lies in the region between the 
curves C_ and C + . Enabling thresholds occur whenever the curve Co intersects the 
curve C_, saturation thresholds occur whenever the curve Cq intersects the curve C+, 
see figure 2. 




Figure 2. The kinematically accessible region for a Lorentz invariant particle of mass m, 
and a Lorentz invariant particle of mass 4m, decaying to two identical particles with energy- 
momentum relation E = \/ Eq + E% tanh (p 2 c 2 /E%). For the mass m particle note the pres- 
ence of both enabling and saturation thresholds. For the mass Am particle only the saturation 
threshold survives. 
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It is additionally useful to distinguish lower and upper thresholds. A lower thresh- 
old occurs when, as a function of increasing po, the curve Co enters the kinematically 
accessible region, and an upper threshold occurs when, as a function of increasing po, 
the curve Co leaves the kinematically accessible region. 




Figure 3. The kinematically accessible region for initial particles with a "lattice- like" 
energy-momentum relation E = (hc/a)\sm(pa/h)\, (with two distinct values of the "lattice 
spacing" a), decaying to two identical particles with "tan/i-like" energy-momentum relation 
E = s/El + El tanh (p 2 c 2 /E 2 ). For small lattice spacing note the presence of four thresh- 
olds: In order they are lower enabling, upper saturation, lower saturation, and upper enabling 
thresholds. For larger lattice spacing only the lower enabling and upper enabling thresholds 
survive. 



There are some quite general results for the curves C± that are not too difficult to 
establish. Consider for simplicity a 2-particle final state, or a 2-particle subsystem of a 
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Figure 4. The kinematically accessible region for a particle with a mass m and a polynomial 
energy-momentum relation, decaying to two identical Lorentz invariant particles of mass m. 
Note presence of both lower and upper enabling thresholds, but no saturation thresholds. 



n-body final state. Then consider the (restricted) set of curves 

C e {p) = {E^pop) + E 2 {[1 - e] Po p), Po } Ve E R, (4.17) 

and (now with p ■ p± = 0) the more general set of curves 

C e>5 (p,p±) = {E 1 (ep p + 8p 1 _)+E 2 ([l-e}pop-8p 1 _), Po } \fe,5 eR. (4.18) 

All of these curves lie between C_ and C+, and can be used to quickly sketch out the 
kinematically allowed region. In particular the curve 

{E^pop/2) + E 2 (p p/2),p } (4.19) 
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Figure 5. The kinematically accessible region for a particle of mass m with a complicated but 
monotonic energy-momentum relation, decaying to two identical Lorentz invariant particles 
of mass m. Note presence of both lower and upper enabling thresholds, but no saturation 
thresholds. 



corresponds to sharing momentum equally between the two particles, and so automat- 
ically lies between C_ and C + . Perhaps less obviously the two curves 

^E 1 (p p) + E 2 (0),p ^ and {^(0) + E 2 (p Q p),p } (4.20) 

correspond to putting all available 3-momentum into particle 1 or particle 2 respectively, 
and both these curves automatically lie between C_ and C + . By considering the limits 
e — > ±oo we also see that the two horizontal lines 

{#i(±oop) + E 2 (Toop), Po } (4.21) 

lie in the kinematically allowed region. (These last two curves are most useful when 
the energy-momentum relations saturate at large 3-momentum.) Finally note that the 
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vertical line 



{E 1 (pop) + E 2 (-p p),o'} (4.22) 



also lies entirely within the kinematically allowed region. These observations allow 
one to quickly sketch key features of the kinematically allowed region. Some graphical 
experiments will quickly convince one that in general the kinematically allowed region 
need not be convex, nor need the curves C± necessarily be monotonic. If one is willing 
to make more specific assumptions concerning the energy-momentum relations, only 
then can much more be said about C± and E max / min (p ) . For instance: 

• For the lattice energy-momentum relation of equation (2.3) we have E^pi) < 
he/ a, so provided all final state decay products see the same lattice, we have 
-E'max(po) — nhc/a, and so the curve C + will be nontrivial. 

• For the tanh-type energy-momentum relation of equation (2.5) it is easy to check 
that 



£max(P0) = E V E h + E li- ( 423 ) 

Then C + is a simple horizontal line. C_ is however quite nontrivial, see figure 2. 

Once one abandons isotropy — in particular azimuthal isotropy around the chosen 
direction p in momentum space — then a fuller analysis using the curves C et s(p,p±) 
will be necessary. Formally 



C-(p) = minminminC e .5(p,p_i_), (4.24) 

eSK pj_ ' ' 



C + (p) = maxmaxmaxC £ .i(p,pi). (4-25) 



and 



Typically the curves C±(p) will piecewise consist of segments of some specific curves 
chosen from the C e ^{p,p±)- The key message to extract from the discussion is this: 
Once exact Lorentz invariance is lost the kinematically allowed region can become 
extremely complicated. 

4.3 Lagrange multiplier techniques 

Motivated by the more restricted analysis of Coleman and Glashow [23], which we now 
significantly generalize, further technical progress can best be made by introducing 
Lagrange multipliers and considering extrema (at fixed po) of the function 

n / n \ 

£( Po ; Pl , A) = E i(Pi) ~ X ■ X> " Po • ( 4 -26) 

i=l \i=l / 
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All minima used to determine E min (p ), or maxima used to determine E max (p ), wm 
be extrema of the function £(p ;pi, A) (though not necessarily vice versa). So extrema 
of S(po;pi, A) will provide information concerning thresholds. 
But all extrema of £ satisfy 

^ = A. (4.27) 
dpi 

In view of the specific Hamilton equation v = x = dE/dp, this implies that at any 
extremum Vi = A = v out — all output velocities for the decay products are equal at any 
extremum, so in particular all output velocities are equal for the specific configuration 
of decay product 3-momenta Pi that define E min (p ) and E max (p ). That is, we have 
the very general result that at threshold all final state particles move with the same 
3-velocity. Furthermore for any extremum we also have 

dS(p ;pi,X) 



dp 



^out, (4.28) 



so in particular this will also be true for _E min (p ) an d E majX (p ). That is: For the 
specific configuration of decay product 3-momenta pi that define .E'min(Po) and E nmx (p ) 
we have 

dp - ° ut ' min ~ d Pi ' [ yj 

and 

dE max dEi 

= «out,max = , (4.30) 

dpo dpi 

respectively. Physically this implies that at any threshold (regardless of whether it is an 
enabling threshold or a saturation threshold, or a lower or upper threshold) all decay 
products will be moving at the same physical 3-velocity. This does not necessarily 
imply that the 3-momenta be related in any simple way, in general the 3-momenta 
need not even be collinear. At threshold we can define the incoming 3-velocity as 

^ = v iQ , (4.31) 
dp 

but with the techniques currently at hand there is in general no simple relation between 
v in and u ut- The best we can currently do is this: If we look along a particular direction 
p in 3-momentum space (with p = p p) then: 

• At a lower enabling threshold p • Vin(po) ^ P ' v out(Po)- 

• At an upper enabling threshold p • Uin(Po) < P • v out(Po)- 
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• At a lower saturation threshold p ■ v in (p ) < p ■ v out (p ). 

• At an upper saturation threshold p ■ v in (p ) > p ■ v out (p ). 

• If p ■ v- m (p ) = p ■ v out (p ) then the curve C (p) touches the kinematically allowed 
region tangentially. One should look at higher derivatives to determine the nature 
of the threshold. If the curve C (p) touches the kinematically allowed region only 
at an isolated point, then we would hesitate to call this any kind of threshold. 
(These isolated points could nevertheless be interesting in their own right.) 

In the case of an isotropic energy-momentum relation this discussion simplifies. The 
momentum p is then parallel (or at worst anti-parallel) to the velocity v and so: 

• At a lower enabling threshold i>i n (po) > ^out(f>o)- 

• At an upper enabling threshold i>i n (po) < v out (po). 

• At a lower saturation threshold i>i n (po) < ^out(po)- 

• At an upper saturation threshold i>i n (po) > ^out^o)- 

• The special case v in (p ) = v out (p ) should be analyzed carefully by looking at 
higher derivatives. This might correspond to an "isolated point at which the 
decay is allowed" ; we would then hesitate to call this any kind of threshold. 

4.4 Thresholds in terms of energy 

With some additional technical machinery we can rephrase the decay thresholds in 
terms of energy rather than 3-momentum. Some aspects of the analysis are more 
complicated, but we will now be able to deduce (at threshold) that v [n and t; out are 
parallel/anti-parallel. Let us now define 



which is the set of all possible total 3-momenta given the input energy E . Now consider 



This is the set of planes (affine subspaces) of co-dimension 3 in IR 3n consisting of all 
possible output 3-momenta compatible with the specified input energy E . Then 




(4.32) 



^out(^in(^o)). 



(4.33) 



E out (V ont (V iQ (E ))) 



(4.34) 
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is the set of all possible 3- momentum-conserving output energies for input energy E Q . 
This will be a connected interval in R. The decay process is then kinematically allowed 
if and only if 

E e E out (V ou t(V in (E ))). (4.35) 

That is 

^min(Pout(^in(^0))) < E < E max (V out (V in (E ))) , (4.36) 
where by this we mean 

min E min (p ) < E < max E max (p ). (4.37) 

Notice now that we are also extremizing over the 3-momenta po compatible with the 
fixed initial energy E Q . Enabling thresholds will then occur at 

min E min (p ) = E Q , (4.38) 
poev in (E ) 

and saturation thresholds at 

E = max E max (p ). (4.39) 

Either one of these thresholds can be characterized in terms of extrema of the related 
function 



£(E ;pi,p ,\,() = ^Eifrd-X 



i=l 



i=i 



C[E (p )-E ], (4.40) 



where we now introduce two Lagrange multipliers, A and (, and we extremize over 
(Pi,Po, A, C) while keeping E fixed. Extremality with respect to the Pi yields 

Vi = A = « out , (4.41) 

whereas extremality with respect to p yields 

C v in = A. (4.42) 

Since the sign and magnitude of ( is unconstrained, this implies that at threshold the 
input velocity v in is either parallel or anti-parallel to the common v out of all the output 
particles: 

C Vm = Vout- (4.43) 

Note that we have gotten at least this far without assuming either spherical symmetry 
or any form of monotonicity. 
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4.5 Asymmetric thresholds 

A particularly peculiar feature of Lorentz violating thresholds is the potential occur- 
rence of asymmetric thresholds, where two identical decay particles might at threshold 
have unequal 3-momenta while traveling at the same 3-velocity. (This phenomena was 
noted, in a more limited context, in reference [27].) The point is that while x(p) is by 
assumption well defined, the inverse function p(x) may be multi- valued. If this happens 
at threshold then two identical particles in the decay channel will have equal velocities 
but unequal momenta. Of course such multi-valued behaviour implies a multi-valued 
Lagrangian L(x) = p{x) ■ x — E(p(x)), which one may wish to exclude from any fun- 
damental theory on physical grounds. (Such behaviour in an effective field theory is 
not particularly problematic.) To characterize when this can and cannot happen, note 
that local invertability of x(p) requires the Jacobian matrix 

I <"•> 

to be nonsingular. Equivalently the Hessian matrix 

dp dp ^ ^ 

should be nonsingular. 

Global invertability of x(p) requires global non-singularity of the Hessian matrix. If 
we now add the extremely mild constraint that the Hessian matrix be positive definite 
at zero momentum (which is required to have any sensible Newtonian or Lorentzian 
limit at low momentum) then global invertability of x(p) requires the Hessian matrix 
to be globally positive definite. But a globally positive definite Hessian matrix implies 
convexity of the energy- momentum relation E(p). 

Thus the existence (or not) of asymmetric thresholds is ultimately related to fail- 
ures (or not) of the convexity of the energy-momentum relation E{p). (For isotropic 
energy-momentum relations, this condition was phrased in terms of a positive curvature 
condition in reference [27].) This is why we can never get asymmetric thresholds in stan- 
dard (non-tachyonic Lorentz invariant) special relativity, and why we do run the risk 
of asymmetric thresholds with (for example) lattice-type, tanh-type, and polynomial 
or rational polynomial energy-momentum relations. Thus the asymmetric threshold 
phenomena encountered by Mattingly, Jacobson, and Liberati in reference [27] is seen 
to have much wider applicability than the situations they considered. 

4.6 Some examples 

As an example of what can happen with asymmetric thresholds, it is quite possible 
for two identical particles to be emitted with almost all the momentum going into one 
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particle, and almost none into the second particle. For instance if one takes two decay 
product particles obeying the tanh-type energy-momentum relation of (2.5), and shares 
the input momentum po in the fractions (| ± e)po, then the final state energy is 

E( Po ,e) = J El + tanh ( [ ^ + j 2 V °^ J + Je 2 + El tanh ( [ ^ "g^ ) . (4.46) 

Depending on the precise ratios between E , E*, and po, this can be minimized at e = 
or near e = 1/2. See figure 6. 




4 - 



-1 -0.5 5 0.5 J 



Figure 6. Energy as a function of e for equation (4.46) holding Eq and E* fixed and for 
five distinct values of pq. Note how the location of the minimum (and hence the threshold) 
shifts from e = (a symmetric threshold) to e ~ ±1/2 (an asymmetric threshold). When the 
minimum occurs at e ~ ±1/2 almost all of the output momentum goes into one of the two 
identical particles, and almost none into the other. 



Another highly nontrivial example, based roughly on equation (3.3), is to take 



E=JEi + c\pl +P 2 y ) + ^ {{Pi - Pi) 3 - Pi}- (4-47) 

This energy momentum relation is carefully chosen to be isotropic at low momentum, 
to have nice behaviour in the x and y directions, and to behave "interestingly" in the 
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z direction. Consider now a particle that moves in the x direction so its original mo- 
mentum is (poj 0, 0). Let it now decay into two identical particles of the type discussed 
above. Imposing 3-momentum conservation, the energy of the final state will be mini- 
mized when the final state particles have 3-momenta (po/2, 0, ±p*). So at threshold the 
initial state and two final state 3-momenta are pointing in three different directions. 
The 3-velocities of the two final state particles will however be equal, and will point 
along the x axis. 

These examples should be viewed as illustrations of the unusual phenomena that 
can occur once strict Lorentz invariance is violated. 

4.7 Monotonicity 

What, if anything, can we say about monotonic energy- momentum relations? (This is 
a common but not universal simplifying assumption.) Consider a set of fixed directions 
Pi and take 

E i (p i )=E i (p i p i ). (4.48) 

For each individual direction pi we can define can define monotonicity in terms of the 
magnitude pi. Monotonicity means 

dE f > 0, (4.49) 
dpi 

but by the chain rule this implies 

^■A>0. (4.50) 
dpi 

That is 

Vi(pi)-pi>0. (4.51) 
But then, by our previous arguments 



dE mi 



dp 



nm/max _ ^ = ^ _ ^ = ^ _ / ^ ^ \ = ^ ^ ■ p { = ^ Vi ■ p { > . (4.52) 



i=l / i=\ 



That is, as long as the individual Ei(pi pi) are monotonic functions of the p i: then 

E m in I max yPo 

Po) is also monotonic as a function of p , and so the curves C± bounding 
the kinematically allowed region will be monotonic. (Monotonicity of the boundary 
curves C± can fail, and quite often will fail, if even one of the final state particles has 
a non-monotonic energy- momentum relation.) 
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4.8 Isotropy 

If all the energy-momentum relations are isotropic (in the preferred (aether) frame) 
then 

Vi oc pf, and v- m oc po- (4.53) 

So all the 3-momenta Pi and po are either parallel or anti-parallel to their corresponding 
3-velocities at threshold. Consequently, in view of the more general results deduced 
above, 

Pi oc p . (4.54) 
That is, all 3-momenta are either parallel or anti-parallel to each other at threshold. 

4.9 Monotonicity plus isotropy 

Only if we assume both isotropy and monotonicity can we deduce that the proportion- 
ality constants in the previous subsection are positive. In this case all 3-momenta and 
3-velocities are parallel at threshold. (This particular theorem was proven by Mattingly, 
Jacobson, and Liberati in [27].) 

5 Scattering thresholds 

Much of the previous discussion of decay thresholds carries over into the discussion 
of scattering thresholds, but there are just enough differences to make some separate 
discussion worthwhile. 

5.1 2-particle collisions 

Consider a 2-particle scattering process of the form 

X A + X B ^X 1 + X 2 + ...X n (5.1) 

involving particles of incoming momenta Pa and p B . Set po = Pa + Pb- We can still 
define both V ou t(Po) and E out (V ou t(Po)), and so construct both -E min (po) and E max (p ). 
The scattering is kinematically allowed if and only if 

E A (p A ) + E B (p B ) e E out (V out (PA + Pb)). (5.2) 
That is, it is kinematically allowed if and only if 

E min (p A + Pb) < E a (pa) + E B (p B ) < E m£bX (p A + Pb)- (5.3) 
Thresholds occur at the boundaries of these regions, that is, at: 

E min (p A + Pb) = E A (p A ) + E B (p B ), (5.4) 
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and at 

Ea{pa) + E B (p B ) = E max (p A + p B ). (5.5) 
To make further progress let us now define 

V in (E A ,E B ) = {po =Pa+Pb- E A {p A ) = E A ,E B (p B ) = E B }, (5.6) 

which is the set of all possible total input 3-momenta given the input energies E A and 
E B . In terms of the notation (4.32) developed for decay processes we can write 

V in (E A , E B ) = V- m {E A ) + V- m {E B ). (5.7) 

Now consider 

Vout(V- m (E A ,E B )). (5.8) 

This is the set of planes (affine subspaces) of co-dimension 3 in IR 3n consisting of all 
possible output 3-momenta compatible with the specified input energies E A and E B . 
Then 

E out (V out (V in (E A ,E B ))) (5.9) 

is the set of all possible 3-momentum-conserving output energies for input energies E A 
and E B . This will be some connected interval in IR. The 2-particle scattering process 
is then kinematically allowed if and only if 

E A + E B e E out (V out (V in {E A , E B ))). (5.10) 

That is 

E min {Vout(V- m {E A , E B ))) <E A + E B < E m .^(V out (V- m {E A , E B ))), (5.11) 
where by this we mean 

min E min (p ) < E A + E B < max E max (p ). (5.12) 

PoeV in (E A ,E B ) PoeV in (E A ,E B ) 

Thresholds will then occur at the edges of the kinematically allowed region. 
Specifically, enabling thresholds will then occur at 



min E min {p ) = E A + E B , (5.13) 

PoeV in (E A ,E B ) 



and saturation thresholds at 



E A + E B = max E max (p ). (5.14) 

Po&VUE A ,E B ) 
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Either one of these thresholds can be characterized in terms of extrema of the related 
function 

E(E A , E B ; pi, p A , p B , X, ( A , Cb), (5.15) 
where we now introduce three Lagrange multipliers, A, (a, and Cb, and set 



E = YE i (p i )-\- 



i=i 



n 



i=l 



(a [E a {Pa) - E A ] - (b \E b {Pb) ~ E B ] . 

(5.16) 



We now hold (E a ,Eb) fixed, and extremize with respect to {p%,Pa,Pb, X, (a, Cb)- 
Extremality with respect to the Pi yields 

Vi = X = v out , (5.17) 

whereas extremality with respect to pa and ps yields 

(a v a = X = Cb v b . (5.18) 

Since the signs (and magnitudes) of (a and Cb are unconstrained, this implies that at 
threshold the two input velocities va and Vb are either parallel or anti-parallel, both 
to each other and to the common v out of all the output particles: 

Ca v a = -«out = Cb v B - (5.19) 

Note that we have gotten at least this far — parallel/anti-parallel input velocities - 
without assuming either spherical symmetry or any form of monotonicity. Without 
additional assumptions we can go no further. 

5.2 Incoming 3-velocities 

In reference [27] Mattingly, Jacobson, and Liberati argue that assuming spherical sym- 
metry and monotonicity of the energy-momentum relations the incoming 3-velocities 
must actually be anti-parallel at (enabling) threshold. Note that in the absence of 
exact Lorentz invariance spherical symmetry can at best only be expected to hold in 
the preferred (aether) fame. 

Certainly if we assume spherical symmetry and monotonicity this result is now 
simple: Assuming spherical symmetry of the individual energy-momentum relations 
Ei(Pi) = Ei{pi) for the decay products, we have E min (p ) = E min (p ). Furthermore 
assuming monotonicity of the individual Ei(pi pi) we have already seen that this implies 
monotonicity of the -E m i n (po p). Then for enabling thresholds the quantity 

min £ min (po) (5.20) 

Po£V ia (E A ,E B ) 
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is minimized when ||po|| i s minimized. Now assuming spherical symmetry for the 
incoming particles, so that Ea/b(pa/b) = Ea/b(pa/b), this occurs when when pa and 
Pb are anti-parallel. Assuming spherical symmetry further implies E(p) = }\\p 2 ) so 
v = dE/dp = f'{\p 2 )p. That is, v and p are either parallel or anti-parallel. But 
monotonicity in turn implies f'(^p 2 ) > 0. That is, the individual v and p are parallel. 
This, in turn, implies the two incoming 3- velocities va and vb are anti-parallel at 
enabling thresholds. 

In contrast, for saturation thresholds we need to consider 



max E max (p ) 



p Q ev in (E A ,E B ) 



(5.21) 



Assuming spherical symmetry and monotonicity for the decay products this quantity 
is now maximized when ||po|| is maximized. But following the argument above, by 
assuming spherical symmetry and monotonicity for the input particles, this in turn 
implies that the incoming 3-momenta Pa and ps are parallel at saturation thresholds, 
which in turn implies that the incoming 3- velocities va and Vb are parallel at saturation 
thresholds. 

Both spherical symmetry and monotonicity are essential to these results. That is: 
Merely deducing that the incoming velocities are parallel/antiparallel at threshold is 
a generic result common to all Hamiltonian-based particle kinematics in a (Lorentz 
violating) homogeneous spacetime. To go further and assert that incoming 3-velocities 
are anti-parallel at enabling thresholds, and parallel at saturation thresholds, requires 
the very much stronger assumptions of spherical symmetry and monotonicity. Our 
results are briefly summarized in tables 1-3. 



Summary of threshold behaviour: Final state particles 





generic 


isotropic 


isotropic+monotonic 


3-velocities 
3-momenta 


equal 
uncorrelated 


equal 
collinear 


equal 
parallel 



Table 1. Behaviour at threshold for the 3-velocities and 3-momenta of outgoing final-state 
particles (compared to each other) under various assumptions. 



6 Conclusions 

As we have seen, abandoning Lorentz invariance carries a very high price. The kinemat- 
ically allowed region, and consequent threshold structure coming from the boundaries 
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Summary of threshold behaviour: Decay — Initial state particle. 





generic 


isotropic 


isotropic+monotonic 


3-velocity 
3-momentum 


collinear 
uncorrelated 


collinear 
collinear 


collinear 
parallel 



Table 2. Behaviour at threshold for the 3-velocity and 3-momentum of the initial decaying 
particle (as compared to the final state decay product particles) under various assumptions. 



Summary of threshold behaviour: Scattering — Initial state particles. 





generic 


isotropic 


isotropic+monotonic 


3-velocities (enabling) 
3-momenta (enabling) 


collinear 
uncorrelated 


collinear 
collinear 


anti-parallel 
anti-parallel 


3-velocities (saturation) 
3-momenta (saturation) 


collinear 
uncorrelated 


collinear 
collinear 


parallel 
parallel 



Table 3. Behaviour at threshold for the 3-velocities and 3-momenta of the two initial state 
(incoming) particles (as compared to each other) under various assumptions. 

of the kinematically allowed region, for both decay processes and 2-particle (elastic 
or inelastic) scattering, is much more complicated than in the Lorentz invariant case. 
There are some limited number of truly general statements that one can make, but 
most of one's intuition has to be reassessed on a case by case basis. We have tried to 
carefully delineate exactly which assumptions are central to which results, concentrat- 
ing on those results that depend only on the existence of a homogeneous spacetime, and 
adding extra assumptions only when essential to obtaining specific specialized results. 

The resulting framework is useful both in (an extremely wide class of) Lorentz 
violating extensions of the standard model of particle physics, and is also potentially 
of interest in quasi-particle settings where violations of Lorentz invariance (and even 
rotational invariance) are the norm rather than the exception. 

Note that in the spirit of classical particle physics we have taken the Hamiltonian 
framework as being more fundamental, and the Lagrangian framework as derivative. 
As a side effect, the Lagrangians that typically arise for non-Lorentz-covariant free- 
particle Hamiltonians are often quite messy and unnatural. If one adopts the view 
that it is the Lagrangian framework that should be viewed as being more fundamental, 
this suggests that modified energy-momentum relations should be most naturally in- 
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terpreted as effective phenomena rather than as fundamental physics. It is ultimately 
for this reason that one commonly focusses on perturbative deviations from Lorentz 
symmetry. A central theme of this article is that such a simplification is not always the 
most useful thing to do, and that there is merit to analyzing thresholds in as general 
a setting as possible. While the analysis of this article has no direct bearing on the 
OPERA-MINOS observations, any serious attempt at phenomenological analysis of 
those observations will need to adopt a theoretical framework along the lines we have 
presented above. 
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